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REACTION-DIFFUSION EQUATIONS 



HANS ENGLER 

Abstract. The fractional reaction diffusion equation dtu + Au = g[u) 
is discussed, where ^ is a fractional differential operator on R of order 
a £ (0,2), the function g vanishes at = and C = 1 a-nd either 
g > on (0, 1) or <; < near — 0. In the case of non-negative g, 
it is shown that solutions with initial support on the positive half axis 
spread into the left half axis with unbounded speed if g{(^) satisfies some 
weak growth condition near = in the case a > 1, or if g is merely 
positive on a sufficiently large interval near ^ = 1 in the case a < 1. 
On the other hand, it shown that solutions spread with finite speed if 
g'{0) < 0. The proofs use comparison arguments and a new family of 
travelling wave solutions for this class of problems. 



1. Introduction 
The scalar reaction-diffusion equation 

(1) dtu{x, t) — dlu{x, t) = g{u{x, t)) 

has been the subject of long study, beginning with the celebrated paper [11]. 
The authors of [11] proposed this equation, with g positive and concave on 
(0, 1) such that g{0) = g{l) = 0, as a model for a population that undergoes 
logistic growth and Brownian diffusion. If u{x,0) = H(x), the Heaviside 
function, and g{u) = u — v?, the equation in fact has an exact probabilistic 
interpretation, see [13] ■ For this problem, is known that solutions approach 
a wave profile ip in the sense that 

u{x + m{t),t) — > ^p{x) {t — > oo) 

where m{t) is the median, u{m{t),t) = ^. It turns out that m{t) = c*t + 
O(logt) for a suitable asymptotic finite wave speed c*. Larger asymptotic 
speeds are only possible if the initial data are supported on M. A more 
general result, given in [2], implies that there is a critical speed c* such that 
for fairly general initial data n(-,0) that are non-negative and supported on 
(0,oo), 

(2) lim limsupii(x, t) = 

i^oo x<~ct 

whenever c > c* and 

(3) lim liminfn(x,t) = 1 

t—>oo x>—ct 
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whenever c < c*. If u denotes a quantity that is to be avoided and whose 
spread is governed by ([T]), a runner may escape from it by running to — oo 
at a speed c > c*, but this quantity will catch up with and engulf her if her 
speed is c < c*. 

Equation ([T|) was also derived in [Ij to describe antiphase domain coarsen- 
ing in alloys. In this situation, g(0) = g{l) = g{u*) = for some u* G (0, 1), 
and g < on (0, u*), g > on (n*, 1). In this case there exists exactly one 
wave speed c* with associated wave profile. In particular, ([2]) still holds for 
this c*. The first of the two cases (the KPP case) corresponds to "pulled" 
fronts (the state u = is unstable) while the second case (the Allen-Cahn 
case) results in a "pushed" front (the state u = is stable). More on these 
two fundamentally different situations may be found in [8] and the refer- 
ences given there. A vast range of applications leading to related models is 
discussed in [9]. 

The purpose of this note is a study of the fractional reaction-diffusion 
equation 

(4) dtu{x,t) + Au{x,t) = g{u{x,t)) . 

Here A is a pseudo-differential operator with symbol p that is homogeneous 
of degree a € (0,2], such that p{—X) = p{X) and \p{l)\ = 1. Then we can 
write p in the form 

(5) p{X) = e-^^'^9n{\)p\xf 

where p £ M. There will be additional restrictions on the parameter p in 
section 2. For p = we obtain fractional powers of the usual negative one- 
dimensional Laplacian, abbreviated often by (— A)"/^. There are various 
real variable representations of such operators, e.g. as a singular integral 
operator or as a limit of suitable difference operators; see [3] where this 
is explained in more detail. The function g is always assumed to satisfy 
g{0) = g{l) = 0. We are interested in both the KPP-case, i.e. g{Q > for 
< C < 1, and the Allen-Cahn case, i.e. g{C) < for ^ near and g{C) > 
for C near 1. 

This class of equations has recenty been proposed as a model for reaction 
and anomalous diffusion; see [31 El [T21 [131 [13 [2D]. It should be noted that 
the term "anomalous diffusion" is also used for situations in which the first 
order time derivative is replaced by a fractional order derivative. Another 
generalization of ([T]) consists in allowing time delays; see |17j . These further 
generalizations will not be discussed here. 

There is strong evidence that the equation (|4]) does not admit traveling 
wave solutions if < a < 2 and g is positive and concave on (0, 1). Rather, 
numerical results in ^ and [4j suggest that for initial data that are supported 
on the positive half axis and increase there from to 1, the median satisfies 
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m{t) 



e^* for some c > 0. In [6j, the estimates 



(6) 



lim limsupii(x, i) = 0, lim liminf m(j;, t) = 1 

t— >oo 3,^_gct t—^oo x>-e'i* 



are shown to hold for such initial data whenever c > c* > d, where c* = 
g'{0)/a. Thus the asymptotic speed of spread grows exponentially. For the 
case where g < on some interval (0, u*), the results in ^X5\ and pl3] suggest 
on the other hand that there exist wave profile solutions that move with 
constant speed, although no rigorous proofs are given there. 

The main results of this note take the form ([2]) and Q. It is shown that 
for a large class of right hand sides g that are non-negative on (0, 1), the 
estimate ([3]) holds for all speeds c; that is, the speed of spread is unbounded. 
It is not necessary to assume that ^'(O) > 0, and li a < 1, g may even be 
zero near 0. On the other hand, if g'{0) < 0, then it is shown that ([2]) holds 
for some finite c, that is, there is always a bound on the speed of spread. 
These results are stated and proven in section 3. This is done by employing 
comparison arguments together with a set of travelling wave solutions that 
is constructed in section 2 and that may be of independent interest. Some 
basic existence and comparison results for ([4]) are sketched in section 4. 



In this section, it will be shown that cumulative distribution functions 
of stable probability distributions with parameters a, [i lead to traveling 
wave solutions u{x,t) = U{x + ct) of (H]), for suitable functions g. A two 
parameter family will be constructed for each possible choice of a and p, 
one parameter being the speed c. The main contribution of this section is 
the characterization of the nonlinear function g that is required to make the 
equation hold. 

Consider the probability density function fa/s of a stable probability dis- 
tribution with index of stability a € (0, 2), skewness parameter (3 € (—1, 1), 
scale parameter 7 = 1, and location parameter 6 = 0, where Zolotarev's 
parametrization (B) (see [21]) is used for a ^ 1 and form (C) is used if 
a = 1. The characteristic function (Fourier transform) of /q,^ then is 



2. A Class of Travelling Wave Solutions 



(7) 



A 1-^ e 



X\°'u}(X) 



where 




The corresponding cumulative distribution function is denoted by F^^. It 
is known that fa/3 is positive, infinitely differentiable, and unimodal. Also, 
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as X ^ oo, there are asymptotic representations 

(8) l-F«/3(x) ~ 

i>i 

(9) fa^ix) ~ '^aCjapX~^~^°' 

and as X — > — oo 

(10) F^pix) - J^d,„;3(-x)-^° 

(11) /a/3(x) ~ 5^ad,„^(-x)-l-^'" 

Consider now the "free" equation 

(12) dtu{x,t) + Au{x,t) =0 
where the symbol of A is given by ([5|), with p given by 

P{a — 1 + sign{l — a)) (a ^ 1) 
/? (a = l). 

Throughout the rest of the paper, we shall only consider operators A with 

symbol p in ^ for which p is of the form ()13p for some (3 G (—1, 1)- The 
relation (jl3p between (3 and p will always be assumed. Taking the Fourier 
transform then shows that equation (jl2p has the fundamental solution 

{x,t)^W{x,t) = t~^l^f^p{xt-^l^) 

with initial data VF(x,0) = 6q{x), the delta distribution. In particular, 
VF(x, 1) = fapix)- There is also the special solution 

(x,t) ^ V{x,t) = Faf^ixt-^^") 

with initial data l^(x,0) = H{x), the Heaviside function. The equations 
hold in the sense of distributions, and the initial data are attained in this 
sense. 

From now on, let a, (5 be fixed. For fixed c G M and r > we consider 
the function 



(13) p 



(14) Urii) = F^p (^r-V- 

Set Ucrix^t) = Ur{x + ct) , then 

Aucr{x,t) = ((x + ct)r-i/"^ 

^ ((x + ct)r-i/"-i) /^^ ((x + ct)T-V") 



a 

and 



dtU,r{x,t) = CT-'/'^f^p ((X + rf)T-^/") 
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and therefore 

dtUcr{x, t) + AUcr{x, t) = ((x + Ct)T-^'^^'^) + CT-^/"^ X . . . 

Now (x + ct)r~-^/" = F~^{ucT{x^t)) and consequently 

(15) dtUcr{x,t) +Aucr{x,t) = CT~^/ goiucAx ,t)) + —gi{ucr{x,t)) 

ar 

with 

(16) go{C) = fap{F-^p{0), gi{C) = F-^{C)fa(,{Fj{C)). 

Equation is of the form @, with g{C) = cr"i/"5o(C) + -^9i{C)- 

In the case a = 1 and —1 < /? < 1, everything is exphcit. Let k = cos ^ 
and a = sin ^ . Then by results in [21j , 



11 X — a 
r Ix) = — I arctan 

F-^(C) = f7-KCOt(7rC) 

1 



7rK(l + (x-a)2/K2) 



5o(C) = — sin2(7rC) 

5i(C) = — sin2(7rC) - -cos(7rC)sin(7rC) 
TTK vr 

It remains to characterize the functions go, gi in the general case. 

Proposition 2.1. Let < a < 2, —1 < (3 < 1. The functions go, gi have 
the following properties. 

a) go and gi are infinitely differentiahle on (0, 1). 

h) The function go is positive on (0,1). The function gi is negative on 
(0,Ffj/3(0)) and positive on (Fq,^(0),1). The function 

C^cr-i/"go(C) + — 5i(C) 
ar 

is negative on (0,u*) and positive on ('u*,l), where u* = -Fq,/3(— car""'^/'^"''^). 

c) AsQi 0, 5o(C) = 0{e+'/^) and 5o(l - C) = 0(C'+'/"). 

d) AsC i 0, gi(C) = -aC + 0(C^+V-). C T 1, 9i{C) = a{l - C) + 

0((1-C)^+^/"). 

e) The functions go and gi can be represented as 

d r^^0^i^ , 
5o(C) = ^ fU')<^' 

d f^^p^^^ 



giiO = ^ sf'^^{s)ds. 
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Proof. Property a) follows since fa/s and Fa(3 together with its inverse are 
infinitely differentiable. Property b) is obvious. Properties c) and d) follow 
from the asymptotic expansions ([9]) and (llip . Finally e) can be checked by 
differentiation. □ 

Property e) will not be used in what follows. It should be noted that go 
and gi are of class on [0,1], but are not infinitely differentiable at the 
interval endpoints, except if a = 1. Clearly go and gi do not depend on c or 
r. We are therefore free to form fairly arbitrary linear combinations of go 
and gi by choosing c and r. 

The construction provides travelling wave solutions for ^ for a special 
class of functions for which g € C^([0, 1]), ^(0) = g{u*) = g{l) for some 
u* G (0,1), and g'{0) < 0, g'{u*) > 0, ^'(1) < 0. This raises the possibility 
that ^ possesses travelling wave solutions for more general functions g with 
these properties. 

If the same construction is attempted for the case a = 2, it turns out 
that go and gi are merely continuous on [0,1], with derivatives that have 
logarithmic singularities near C = and C = 1. Therefore the arguments 
in the next section cannot be extended to the case a = 2, and indeed the 
results of the next section do not hold in that case. 



3. Results on the Speed of Spread 

This section contains the main results of this paper. As before, the opera- 
tor A has symbol ([5]) with < a < 2 and p satisfying ()13p with —1 < /3 < 1. 
We always assume that n is a solution of (|4]) and that g G ^""^([0, 1],M) with 
g{0) = g{l) = 0. Initial data uq will be assumed to satisfy 

(17) no G C(M,M), < uo{x) < 1, lim uo{x) = 1, supp{uo) C [0,oo) . 

The results in section 4 then imply that (j3]) has a unique mild solution u 
that exists for all x G M, t > 0, and this solution satisfies < u{x, t) < 1 for 
all {x,t). The notation of that section will also be used here. 

We first discuss the case where g > on (0, 1). The main result in this 
case is the following. 

Theorem 3.1. Let u be the solution of dH) with uq satisfying (jl7p . 

a) Let a > 1. Assume that g > on (0, 1) and that for some cq > 0, < 
7 < ^ and all ( G [0, i] 

9iO>coC- 

Then for all c > 

liminf inf u{x,t) = 1. 

t — ^oo x>—ct 

h) Let a = 1. Assume that g > on (0, 1). Then for all c > 

liminf inf u{x,t) = 1. 

t — ^oo x>—ct 
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b) Let a < 1. Assume that g > on (0, 1) and g{C) > for ( E [2^' -'^)- 
Then for all c > 

liminf inf u{x,t) = 1. 

t^co x>—ct 

The result shows that the speed of spread is unbounded (that is, ([3]) holds 
for all c > 0), and it exhibits different mechanisms for this phenomenon. 
Recall that in the interpretation of [11] , the function g is responsible for the 
growth of a substance whose density is given by n, while A describes the 
spread of this substance. If a G (0,2), the substance spreads with a jump 
process, not with Brownian diffusion, and jumps of magnitude h occur with 
a probability that is 0{h~°') for large h. If a > 1, the mean jump distance 
is still finite. In this case, the growth rate g{u) at small densities (small 
u) is responsible for the unbounded speed of spread. If a is close to 1, 
this growth can be very weak {g{C) ~ O with large 7), yet the speed of 
spread is still unbounded. If on the other hand a < 1, i.e. jump sizes 
have unbounded mean, the growth rate for small densities does not matter 
any more for the speed of spread to be unbounded; in fact there may be 
no growth at all for small densities {g{u) = for small u), and yet the 
speed of spread is unbounded. In this case, the unbounded speed of spread 
results from growth that occurs solely for large densities {g{C) > only for 
C > ^^)- The substance is transported towards —00 due large (a < 1) 
negative jumps, resulting in an unbounded speed of spread. It is known 
that in this case, is the fraction of negative jumps. If this fraction is 
large, then it is sufficient that growth occurs only for densities close to the 
maximal value, i.e. g{u) > on [^-^, 1) already implies that the speed of 
spread is unbounded. The case a = 1 is intermediate: Any growth for small 
densities {g{C) > for C > 0) results in an unbounded speed of spread. 

In the case a > 1, it would be interesting to know if the speed of spread 
is still unbounded if 7 > or if a finite speed of spread occurs holds 
for large c) if 7 becomes sufficiently large, i.e. if growth is extremely weak 
for small densities u. In the case a < 1, it would be interesting to know if a 
finite speed of spread is possible at all ii g > and g is not identically equal 
to 0. 

The main result in the case where g is negative near C = is the following. 

Theorem 3.2. Let u be the solution of ([Ij) with initial data uq satisfying 
(I17p . Assume that g'{0) < 0. Then there exists c > such that 

limsup sup u{x,t) = 0. 

t— ►00 x<—ct 

The result shows that negative proportional growth at small densities 
(g'iO) < 0) always limits the speed of spread of a substance whose growth 
and spread are governed by ([4]), even for processes whose jump sizes tend 
to be very large (a < 1). I am not aware of an interpretation of ([4]) in the 
context of material science, similar to the use of ([T]) in [IJ. 
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The proofs will be given below. The main tools in the proofs are the 
comparison arguments given in the next section, together with the following 
auxiliary result. 

Lemma 3.1. Let g G C^([0, 1],]R) and let g^, gi be defined as in (I16p . 
depending on a G (0,2), /? G (—1,1). 

a) Let a € (1,2). Suppose that g{() > for all ( e (0, 1] and that there 
exist Co > and 7 < such that g{C) > cqC^ for all ( G [0, ^]. Then 
given any c > there exists r > such that for all C G [0,1] 

5(C)>cr-i/"go(C) + («r)-V(C). 

b) Let a = 1. Suppose that g > on (0, 1]. Then given any c > there 
exists r > such that for all ( € [0, 1] 

9{0>cT-'9o{C) + T-'gi{C)- 

c) Let a G (0,1). Suppose that g > on [0,1] and g{Q > for all 
C € [i^,l). Then given any c > there exists r > such that for all 
C G [0, 1] 

<7(C)>cr-i/"5o(C) + («r)-V(C). 

d) Suppose g'{0) < and g{C) = for C € [1 — e, 1] for some e. Then 
there exist c € M and r > such that for all C E [0,1] 

<7(C)<cr-i/"go(C) + (ar)-V(C). 

Proof. Consider first statement a). Let a > 1, —1 < (3 < 1 and let M > 
be large enough such that for some ci, C2 > and all x < — M 

This is possible by (jlip . Let c > be given, then we may increase M further 
such that also 

, ^r- a — 1 _2_ 

CqcIM"^ > C2 c"-i . 

a 

where r = -^^^ — 7 > 0. We omit the subscript a/5 in the formulae involving 
and f^p from now on. Now set 5 = F{—M). Then for < C = F{x) < 5, 
i.e. X < -M, and for ah r > (M/c)"''^""^) 

CT-'/-go{Fix)) + {aT)-'g^{Fix)) = (^cr-^/^ + {arr'x) f{x) 

< ^c^|x|i/(i-")/(x) 
a 

a 

= C'^-^C2\x\ ' 

a 

where a standard calculus argument has been used to see that the expression 
(cT^^/° + (ar)^-'^x) is maximal for r = {\x\/ c)°'^'^°^'~'^\ We estimate further. 
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using the choice of M 

< CoC^lxr'lxr"'/^"-^) = Co (cilxD^ 

< co{F{x)r <g{F{x)). 
Therefore, for all r > (M/c)"/^""^^ and all C < = F{-M), 

cr-i/"5o(C)) + (ar)-V(C)<5(C). 

Since gf > on [5, 1] by assumption, this inequality can be achieved also on 
[5, 1] by increasing r even further. This proves part a). 

The proof of part b) is straight forward: Given c > 0, note that C(7o(C) + 

51 (C) < on the interval [0, F{-c)]. Then ^^oiC) + giiC) ^ ^^^^ ^^^^ ^ 

T 

is sufficiently large. 

To prove part c), let again c > be given. Let u* = inf{n € [0, 1] \g{u) > 
0}. Then u* < = F{0). Pick r large enough such that F (— car^~^/") > 
u*. This is possible since a < 1. Then on [0, u*], 

cr-i/"5o(C)) + (ar)-i5i(C)<5(C) 

since the left hand side is non-positive there by Proposition [2Tj By increas- 
ing T further, we can obtain this estimate also for € using again 
that g is assumed to be positive on [^-^, 1]. 
To prove part d), note first that 

cr-i/"5o(C)) + (ar)-i5i(C)>^(C) 

on an interval [0,6] as soon as + ^'(0) > 0, i.e. for sufficiently small 
r. Increasing c sufficiently and noting that g = near C = 1 extends this 
inequality to the entire interval [0, 1]. □ 

Proof of Theorem 3.1. The proof uses the same argument for all three parts, 
so we give details only in part a). Let e > 0. We replace u with n = (1 + e)u 
and g with g, where g{C) = (1 + e)5((l + e)~^C)- Then 

dtu + Au = g{u) . 

Then g{() > cqC^ for C £ [0, possibly with a changed cq, and additionally 
g > on (0, 1]. Let c > be given, then there exists r > such that 

(18) m > (c + l)r-V"g^(^) + (arr'g.iC) . 

for all C G [0, 1] by Lemma [3.1l By extending go and gi to be zero on [1, l-|-e], 
this inequality is true on [0, 1 + e]. Now find a constant d such that vq{x) = 
H{x — d) < u{x, 0) for all x. This is possible since lima;^oo uix, 0) = 1 -|- e. 
By Proposition I4.3t we see that 

u{x,t) > F (^{x - d)t-^/"^ 
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for all X G M, t > 0. This is in particular true for t = t. Now use Proposition 
lO and ([18]) to infer that 

u{x, t)>F{{x-d+{c+ l)t)r-i/") 

for all x E M, t > r. Therefore for t>T and x > —ct, 

u{x,t) > F(^(x-d+(c+l)t)T-i/°) 

> F((-ct-d+(c+l)t)r-i/"' 



As t ^ oo, the right hand side goes to 1. Rewriting this in terms of u, we 
see that 

liminf inf u{x,t) > {1 + e)~^ . 
t— >oo x>—ct 

Since e was arbitrary, the desired result follows. 

In case of part b), the same argument can be used without changes, 
appealing to part b) of Lemma 13. 1[ 

In case of part c), we have to restrict e such that ^ > on [''-^, 1], that is, 
g > on [-^Y+t)^ "^^^ ^^^^ °^ proof is again unchanged, using part c) of 
Lemma 13.11 □ 

Proof of Theorem 3.2. We replace u with u = and g with g, where g{(^) = 

y{2C) for < C < ^ and ^(C) = for C £ !]• Then ^'(0) = ^'(0) < 
and 

dtu + Au = g{u) . 
By Lemma |3.H part d), there exist c > and r > such that 

5(C)<(c-l)r-i/"5o(C) + («r)-V(C). 
Since lim^j^oo u{x, 0) = | and u{x, 0) = for x < 0, we can find d > such 
that F{x + d)T-^/'^) > 0) for all a; G M. Using Proposition 14. H one 
sees that 

F[{x + d+{c- l)t)T-^/") > u{x, t) 
for all X € M, t > 0. Therefore for i > and x < —ct, 

u{x,t) < F('(x + d+ (c- l)t)T-i/"' 



< F(^(-ct + d+(c-l)t)r-i/° 

The right hand side tends to as t ^ oo. In terms of u, this implies 

limsup sup u{x,t) = 0. 

t— >oo x<—ct 

This concludes the proof. □ 
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4. Facts About Fractional Reaction-Diffusion Equations 

In this section we summarize some basic theory about (j3|) that is needed 
in this note. A broader and deeper discussion may be found in [3]. 
We work in the Banach space 

Ciim = {w £ C(M) I hm w{x) and hm w{x) exist} 

X — ^oo X — ^— oo 

equipped with the supremum norm || • ||. Let A be the pseudodifferential 
operator with symbol ([5]) and parameters a, p. As always, let p, (3 be related 
by (|13p and — 1 < /? < 1. Solutions of the free equation (jl2p with initial 
data u(-,0) = (/9 € Cum then can be written in terms of the fundamental 

solution {x,t) ^ t'^^"faf3 (It^)) namely 

(19) t)=J^ t-'/y^p cp{y) dy 

where fap is a stable probability density function. 
For fixed a and p and ip E Cum-, define 

S{t)ip{x) = t) 

where u is given by ([T9|) . This is a positive Co semigroup on Cum and a Feller 
semigroup on the subspace of functions in Cum that vanish at ±00. If ■0 is a 
continuous function from [0, T] to Cu^^ then solutions of the inhomogeneous 
equation 

(20) dtu{x,t) + Au{x,t) = '4){x,t), u(-,0)=99 
can be written with the variation-of-constants formula 

(21) u{-,t) = S{t)^+ f S{t- s)ip{-,s)ds. 

Jo 

A continuous curve u : [0, T] — > Cum that satisfies (I2ip is commonly called 
a mild solution of (j20p . Next let g : [0,oo) x M — M be locally Lipschitz 
continuous in both variables and let ip S Cum- Then the equation dtu{x, t) + 
Au{x^t) = g{t,u(x,t)) (for which (|4]) is a special case) has a unique mild 
solution u £ C {[0,T),Ciim), where < T < 00 is maximal. Either T = 00, 
or ||u(-,i)|| — > 00 as t t T. The solution can be obtained as the locally in 
time uniform limit of the iteration scheme 

Un+i{-,t) = S{t)Lp + S{t - s)g{s,Un{-,s))ds (n = 0,l,...) 
Jo 

with uq arbitrary, e.g. uo{-,t) = S{t)ip. It is possible to set up a more general 
solution theory, but this is not needed for the purposes of this paper. 

Solutions of (jl]) satisfy comparison theorems. Results of this type are true 
for all Feller semigroup. A systematic study of such semigroups and their 
generators was carried out in [5] , following the seminal work on this topic in 
|18j . For the sake of completeness, a comparison result is stated here, and 
its proof is sketched. 
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Proposition 4.1. Let u, v € C {[0,T],Ciim) be mild solutions of the equa- 
tions 

dtu + Au = g{u), dtv + Av = h{v) 
where 5', /i : M — >■ M are locally Lipschitz continuous. If 

5(0 < MO vcgm 

and 

Ui;0)<vi;0) 

then 

u{x,t) <v{x,t) V(x,t) G M X [0,r] . 

Proof Let M = max[o,T] {\H-,t)\\ + \\v{-,t)\\ + 1). Let A > \g'{C)\ + h'{C)\ 
for all |(^| < M. Without loss of generality we may assume that g and h are 
constant outside [— M, M]. Set 

U{x, t) = e^^u{x, t), V{x, t) = e^^v{x, t) 

and observe that U and V satisfy 

dtU + AU = g{t,U) 
dtV + AV = h{t,V) 

with g{t, = XC+e^^g (e~'^*C) and h{t, Q defined similarly. Clearly, g{t, Q < 
h{t, C) for all (. The function g is non-decreasing in its second argument, 
since for almost all ( 

WO = A + 5'(e-^^)>0. 
Consider the iteration scheme 

Un+i{;t) = S{t)u{;0)+ [ S{t-s)gis,Un{;s))ds 

Jo 

and similarly for and h. The scheme for the Un converges to the limit U, 
and the scheme for the Vn converges to the limit V. 

We now employ a standard induction argument to show that Un <Vn on 
M X [0, T] for all n. This implies that U < V and therefore also u < v on 
M X [0,r]. Let Uoi-,t) = S{t)u{-,0) and Vo{-,t) = S{t)v{-,0), then Uq < Vb 
on R X [0,T] since S is a positive semigroup and u(-,0) < f(-,0). Suppose 
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[/„ < K on M X [0,T], then 

Un+i{;t) = S{t)u{;0)+ [ Sit-s)g{s,Un{;s))ds 

Jo 

< S{t)v{;0)+ [ S{t-s)g{s,Un{;s))ds 

Jo 

< S{t)v{;0)+ [ S{t-s)g{s,Vn{;s))ds 



< S{t)v{-,0)+ S{t-s)his,Vni-,s))ds 
Jo 

= Vn+l(,;t) 

which completes the induction step. This proves the proposition. □ 

Corollary 4.2. Consider a mild solution u € C {[0,T),Ciim) of @ and 
assume that g is locally Lipschitz continuous. If 3(7) > for some 7 and 
u{-,0) > 7, then u{-,t) > 7 for all t. If 3(7) > > g{5) for some j < 6 
and 7 < n(-,0) < 6, then 7 < u{-,t) < 6 for all t, and the solution can be 
continued to M x [0, 00). 

The proof consists in observing that the constant functions v{x,t) =7 and 
w{x, t) = 6 solve with right hand sides and therefore must be pointwise 
bounds for the solution, by Proposition l4.1[ If the solution remains bounded 
between two constants, then its supremum norm remains bounded and it 
can be continued to R x [0, cxd). 

Also required is a comparison result for solutions whose initial data are 
step functions. Since such initial data are not in Cum, a separate argument 
is required. 

Proposition 4.3. Let u G C {[Q,T],Cum) he a mild solution of ([1]), with 
locally Lipschitz continuous g. Assume that 

N 

u{x, 0) > ^0(2;) = oo + ajH{x — Cj) Vx € M 

i=i 

where aj G M, ci < C2 • • • < cn, and H is the Heaviside function. Let 
7 = min]Rfo(3;) and 6 = max]8Wo(x). Assume also that g > on [fjS]. 
Then 

^ / _ \ 

(22) u{x,t)>ao + ^ajF^J^^-^\ VxGM,0<t<r 

i=i 

where F^p is the cumulative distribution function of the associated stable 
distribution. 
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Proof. We know that u{x, t) > 'j and thus may assume that g{C) > also 
for C < 7- For arbitrary ea > 0, we set 

1 r 

Vea{x) = — I Vq{x — z)dz — 6 . 

o" Jo 

Then is piecewise hnear and constant outside the interval [ci , cat + a] ; in 
particular, feo- ^ C^jm- Solving (jl2p with initial data ^ gives the solution 




Given e > 0, it is possible to find o" > such that v^^j < uq{x) on M, since uq 
is uniformly continuous. Clearly, 7— e < V^a < 5- We may therefore view V^a 
as a solution of with a right hand side h that satisfies = < g{zeta) 
for C ^ 1^ S'lid < 5(C) ^-Iso for > 5). Then by Proposition 14.11 

u{x,t) >V,^^{x,t) VxGM, o<t<r. 

Send 5 to 0, then since Fa/3 is uniformly continuous, (j22|) is obtained with 
qq replaced by oq — e on the right hand side. Now send e to and (^Z2\i 
follows. □ 
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